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Framework  for  Control  of 
Dynamic  Ice  Breakup 
by  River  Regulation 

MICHAEL  G.  FERRICK  AND  NATHAN  D.  MULHERIN 


Introduction 

The  breakup  of  a  strong  and  intact  river  ice 
cover  at  a  high  discharge  is  a  dynamic  process 
with  the  potential  to  cause  significant  damage. 
For  example,  Ferrick  et  al.  (1988)  concluded 
from  a  study  of  historical  data  that  the  highest 
water  levels  and  greatest  potential  for  damage 
to,  or  loss  of,  a  historic  covered  bridge  over  the 
Connecticut  River  at  Windsor,  Vermont,  is 
present  during  a  dynamic  ice  breakup.  When 
the  ice  movement  occurs  on  the  Connecticut 
River,  the  resulting  channel  blockage  and  in¬ 
creased  hydraulic  roughness,  combined  with  a 
high  discharge,  can  cause  extreme  water  levels 
that  threaten  the  bridge  and  cause  flood  dam¬ 
age  in  Windsor.  However,  the  data  indicate 
that  the  occurrence  of  large  breakup  events  is 
predictable,  with  10  of  12  recent  large  events 
occurring  in  a  2-1 /2-week  period  in  March, 
each  following  and  in  response  to  a  significant 
rainfall. 

Current  methods  used  to  mitigate  the  effects 
of  dynamic  ice  breakup  include  blasting  of  ice 
jams  that  have  already  formed,  channel  excava¬ 
tions  to  prevent  jam  formation,  construction  of 
dams  or  ice  retention  structures  to  control  the 
movement  of  the  ice,  and  thermal  discharges  to 
melt  ice  prior  to  breakup.  The  thermal  discharge 
method  requires  a  source  of  warm  water  that 
may  not  be  available,  and  the  other  methods 
may  not  be  effective  or  environmentally  accepta¬ 
ble.  However,  in  regulated  rivers  the  capability 
to  manipulate  the  flow  may  be  used  to  control 
both  the  ice  formation  and  breakup  processes. 
Donchenko  (1978)  observed  that  intensive  defor¬ 
mation  or  breakup  of  an  ice  cover  downstream 


of  hydropower  dams  occurs  in  response  to  a  rap¬ 
id  rise  in  river  stage.  He  suggested  that  flow  re¬ 
leases  could  be  regulated  to  minimize  ice  break¬ 
up  and  jamming  during  the  period  of  ice  cover 
formation.  The  development  of  generally  appli¬ 
cable  ice  management  strategies  based  on  river 
regulation  requires  a  quantitative  theory  of  dy¬ 
namic  ice  breakup. 

The  primary  contribution  of  this  work  is  a 
quantitative  description  and  subsequent  demon¬ 
stration  of  the  fundamental  relationship  between 
unsteady  flow  and  dynamic  ice  breakup.  At 
breakup,  the  hydraulic  forces  on  the  ice  cover 
vary  widely  in  response  to  unsteady  flow,  and 
this  variability  cannot  be  neglected  without  care¬ 
ful  consideration.  The  peak  forces  on  the  ice  cov¬ 
er  do  not  necessarily  coincide  with  either  the 
peak  stage  or  the  peak  discharge.  In  this  report 
we  describe  and  classify  the  range  of  ice  breakup 
behavior  as  completely  as  possible,  consistent 
with  our  present  understanding.  Building  on 
this  description,  we  refine  the  theory  and  nu¬ 
merical  model  presented  by  Ferrick  et  al.  (1986b) 
by  developing  a  force  balance  for  a  common  dy¬ 
namic  breakup  behavior.  We  evaluate  the  effect 
on  the  ice  cover  of  these  forces  with  empirical 
breakup  resistance  criteria.  These  criteria  are  de¬ 
veloped  for  a  case  study  of  the  Connecticut  Riv¬ 
er,  using  the  model  together  with  field  data  ob¬ 
tained  during  a  dynamic  breakup,  and  are  com¬ 
pared  with  published  values.  The  force  balance 
can  then  be  applied  to  analyze  the  collapse  of  a 
river  ice  cover,  including  relatively  thin  ice  that 
occurs  during  the  ice  formation  process  and  the 
formation  and  release  of  ice  jams.  An  application 
of  the  completed  model  demonstrates  the  in- 
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sights  for  interpretation  of  observations  that  fol¬ 
low  from  the  theory  presented,  and  the  intuitive 
nature  of  these  results. 

This  framework  for  understanding  river  ice 
processes  provides  the  option  for  ice  manage¬ 
ment  by  flow  regulation.  We  focus  on  the  poten¬ 
tial  for  control  of  ice  breakup  on  regulated  riv¬ 
ers.  Model  sensitivity  studies  indicate  that  the 
ice  cover  response  to  a  controlled  release  can 
vary  greatly  with  ice  thickness  and  breakup  re¬ 
sistance. 

Characteristics  of  ice  breakup 

River  waves  are  long-period,  shallow  water 
waves  that  are  a  consequence  of  unsteady  flow 
(Ferrick  1985).  Flow  releases  at  hydroelectric 
dams  typically  cause  abrupt  river  waves  that  can 
delay  or  prevent  the  formation  of  an  ice  cover 
and  cause  ice  jams,  indicating  an  important  role 
of  these  waves  in  ice  behavior.  The  importance 
of  river  waves  in  the  ice  breakup  process  has 
been  noted  by  several  authors  (Beltaos  and 
Krishnappan  1982,  Billfalk  1982,  Doyle  and  An¬ 
dres  1979,  Henderson  and  Gerard  1981,  Prowse 
et  al.  1986,  and  Wong  et  al.  1985).  It  is  widely  rec¬ 
ognized  that  ice  jam  formation  and  release  ini¬ 
tiates  river  waves,  and  that  waves  associated 
with  sudden  jam  release  cause  ice  breakup.  The 
hydraulic  forces  on  a  river  ice  cover  are  related 
to  the  flow  energy  gradient,  a  parameter  that  can 
increase  substantially  on  the  front  of  a  river 
wave.  River  ice  breakup  occurs  when  the  forces 
on  the  cover  exceed  the  resistance  provided  by 
the  ice  strength  and  points  of  support.  With  an 
ice  sheet  in  place,  the  frictional  resistance  to  the 
flow  of  a  river  is  increased,  affecting  both  the 
steady-state  stage-discharge  relationship  and  the 
unsteady  flow  dynamics.  River  waves  and  dy¬ 
namic  ice  breakup  are  intrinsically  related  be¬ 
cause  waves  form  during  a  dynamic  breakup 
from  the  release  of  water  in  channel  storage  with 
the  rapid  decrease  in  flow  resistance  as  the  ice 
breaks  up. 

Breakup  may  occur  at  any  point  on  the  spec¬ 
trum  between  small  forces  that  exceed  a  greatly 
diminished  ice  cover  resistance,  characteristic  of 
a  thermal  breakup,  and  very  large  forces  that 
overcome  the  resistance  of  thick  and  competent 
ice,  termed  a  dynamic  breakup.  The  dynamic  or 
thermal  character  of  ice  breakup  on  a  river  will 
typically  span  the  spectrum  over  a  period  of 
years.  Zachrisson  (1988)  reported  widely  differ¬ 
ing  breakup  behavior  in  different  years  of  the 
unregulated  River  Tornealven  at  the  border  be¬ 


tween  Finland  and  Sweden.  High  water  levels 
occur  when  the  ice  is  thick  and  breakup  is  pre¬ 
ceded  by  rapid  river  stage  increases,  characteris¬ 
tics  of  dynamic  ice  breakup.  The  only  condition 
common  to  mild  breakups  on  this  river  is  a  slow 
increase  in  river  stage  before  breakup,  a  charac¬ 
teristic  of  thermal  events.  Prowse  et  al.  (1988) 
performed  pre-breakup  tests  of  in-situ  ice 
strength  on  the  lower  Liard  River  in  northern 
Canada.  They  reported  a  linear  decrease  in  ice 
strength  of  50%  over  an  18-day  period,  followed 
by  a  relatively  mild  breakup.  The  extent  of  dete¬ 
rioration  of  ice  strength,  combined  with  the  mag¬ 
nitude  and  rate  of  change  in  hydraulic  condi¬ 
tions,  determines  the  character  of  breakup.  The 
development  of  abrupt,  high-amplitude  river 
waves  during  breakup  requires  significant  run¬ 
off  and  competent  ice  cover. 

The  characteristics  of  a  dynamic  ice  breakup 
depend  on  the  mode  of  failure  of  the  ice  cover. 
Ferrick  et  al.  (1986b)  referred  to  support- 
dominated  and  strength-dominated  breakup  be¬ 
haviors  as  of  high  and  low  energy,  respectively. 
A  failure  at  the  supports  of  an  ice  cover  produc¬ 
es  a  sudden  bank-to-bank  release  of  the  ice.  This 
support-dominated  breakup  travels  rapidly 
downstream  at  a  speed  greater  than  the  flow  ve¬ 
locity.  The  breaking  front  is  the  boundary  be¬ 
tween  the  intact,  stationary  ice  cover  down¬ 
stream  and  moving  ice  plates  upstream  (Fig.  1). 
Fractures  in  the  ice  cover  at  preexisting  cracks 
appear  with  the  initial  motion,  resulting  in  a 
change  from  a  continuous  cover  to  ice  plates. 
The  high  speed  of  the  breaking  front  relative  to 
that  of  the  ice  plates  prevents  ice  participation  in 
the  breakup  downstream.  The  sizes  of  the  mov¬ 
ing  plates  are  reduced  with  time  due  to  colli¬ 
sions.  The  rubble  front  separates  the  ice  plates 
from  the  brash  ice,  and  represents  a  region  of  ice 
convergence  that  can  develop  a  significant  thick¬ 
ness.  Finally,  open  water  predominates  behind 
the  ice  run. 

We  have  observed  an  arrest  of  the  motion  of  a 
support-dominated  breaking  front  at  locations 
with  reduced  energy  gradient,  and  a  transition 
to  a  strength-dominated  breakup.  The  rubble 
front  and  the  breaking  front  are  coincident  in  a 
strength-dominated  dynamic  breakup.  This 
breaking  front  initially  extends  across  only  the 
high  velocity  portion  of  the  channel,  and  the  ice 
in  more  sheltered  locations  moves  somewhat  lat¬ 
er.  Brash  ice  from  upstream  interacts  with  the 
leading  edge  and  is  transported  under  the  intact 
cover.  At  the  front,  relatively  small  pieces  of  ice 
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Rubble  Front 


B'eak  ng  Front 


Figure  1.  Overview  of  a  support-dominated  dynamic  ice  breakup.  The 
breaking  front  travels  downstream  at  a  speed  of  several  meters  pet- 
second  and  is  distinct  from  the  ice  rubble  front.  This  depicts  a  single 
breaking  front  in  a  progressive  breakup. 


break  in  bending  from  the  intact  sheet.  The 
progress  of  a  breakup  that  requires  the  ice  to  fail 
is  slower  than  that  of  a  breakup  characterized  by 
support  failure.  Strength-dominated  breakup 
can  occur  at  significantly  lower  hydraulic  forces 
than  support-dominated  breakup,  and  the  break¬ 
up  celerity  is  directly  related  to  the  magnitude  of 
these  forces.  The  breaking  front  speed  is  the  ve¬ 
locity  of  the  brash,  indicating  the  importance  of 
the  interaction.  If  the  speed  of  breakup  is  signifi¬ 
cantly  lower  than  the  river  wave  speed,  a  gradu¬ 
al  separation  develops  between  the  wave  front 
and  the  breaking  front,  and  the  moderate  rate  of 
water  release  from  channel  storage  may  not  be 
sufficient  to  prevent  wave  attenuation.  Both  of 
these  factors  cause  the  hydraulic  forces  on  the 
ice  to  diminish  with  time,  eventually  leading  to  a 
complete  stall  of  the  breaking  front  and  the  for¬ 
mation  of  an  ice  jam.  The  release  of  this  jam  ini¬ 
tiates  the  next  surge-stall  sequence,  a  behavior 
that  is  characteristic  of  dynamic  ice  breakup. 

The  concept  of  celerity  matching  between  the 
breaking  front  and  the  river  wave  is  needed  for 
understanding  the  sequencing  of  dynamic  break¬ 
up  behaviors.  When  these  celerities  are  approxi¬ 
mately  equal,  the  hydraulic  forces  driving  the 
breakup  are  maintained  at  a  high  level.  A  slow¬ 
ing  of  the  breaking  front  relative  to  the  wave, 
due  to  increased  breakup  resistance  at  some  lo¬ 
cation,  typically  corresponds  to  a  diminished 


force  environment,  and  further  slowing  or  stall¬ 
ing  of  the  breakup  is  possible.  The  outcome  de¬ 
pends  on  the  duration  of  the  slowing  relative  to 
the  wave  celerity,  the  displacement  of  the  break¬ 
up  from  the  position  on  the  wave  of  the  highest 
forces,  and  the  magnitude  of  the  force  reduction. 
Sustained  and  rapidly  moving  strength-domi¬ 
nated  breakups  are  possible  in  rivers  that  are 
steep  and  hydraulically  rough  because  of  the 
dominance  of  bulk  waves  (Ferrick  1985).  The 
reason  is  that  the  bulk  wave  celerity  is  approxi¬ 
mately  equal  to  the  flow  velocity,  and  if  a 
strength  breakup  has  this  same  celerity,  the 
matching  can  cause  sustained  movement  over 
large  distances.  In  mildly  sloped  rivers,  celerity 
matching  at  the  flow  velocity  does  not  occur  and 
strength  breakups  at  this  speed  are  not  stable. 

We  use  the  term  progressive  to  indicate  a  dy¬ 
namic  breakup  advance  where  the  front  of  the 
breakup  moves  downstream  coincident  with  the 
front  of  a  single  dominant  river  wave.  A  pro¬ 
gressive  breakup  can  exhibit  both  strength-  and 
support-dominated  failure,  and  surge-stall  se¬ 
quences.  Flood  water  levels,  scour  of  river 
banks,  and  damage  to  structures  near  the  river 
are  common  with  progressive  dynamic  break¬ 
ups.  In  contrast,  a  simultaneous  breakup  occurs  at 
several  locations  concurrently  and  is  associated 
with  weakened  ice  conditions  or  with  hydraulic 
forces  that  are  marginally  adequate  to  produce  a 
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Figure  2.  Stage  record  of  a  progressive  ice  breakup  of  the  unregu¬ 
lated  White  River  on  2  February  1988.  This  wave  developed  as  a 
result  of  ice  breakup  upstream.  The  large  ice  forces  associated 
with  the  wave  front  caused  the  breakup  to  move  downstream  with 
the  wave. 


breakup.  Simultaneous  breakups  can  be  classi¬ 
fied  as  thermal,  dynamic,  or  somewhere  in  be¬ 
tween,  depending  on  the  strength  and  thickness 
of  the  ice  cover,  and  the  amplitude  and  steep¬ 
ness  of  the  waves  associated  with  the  breaking 
fronts. 

A  progressive  breakup  can  be  initiated  in  un¬ 
controlled  rivers  by  a  simultaneous  breakup  in 
the  upland  portion  of  a  basin.  Sharp-fronted  riv¬ 
er  waves  that  develop  in  unregulated  rivers  can 
cause  a  progressive  dynamic  breakup  over  great 
distances,  with  few  stalls  of  relatively  short  du¬ 
ration.  Other  important  characteristics  of  these 
waves  are  significant  stage  increase,  short  dura¬ 
tion  of  the  high  water  levels,  and  a  large  contrast 
in  discharge  between  the  initial  and  peak  flow 
conditions.  Figure  2  presents  stage  data  from  a 
progressive  breakup  on  2  February  1988  of  the 
unregulated  White  River  in  Vermont.  Ice  break¬ 
up  occurred  at  this  gaging  station  immediately 
following  the  arrival  of  the  wave.  This  record 


does  not  exhibit  abrupt  rises  and  falls  character¬ 
istic  of  ice  jamming.  The  sharp-fronted,  high- 
amplitude  and  short-duration  features  of  this 
river  wave  are  all  characteristic  of  a  progressive 
ice  breakup. 

The  characteristics  of  dynamic  ice  breakup  on 
regulated  rivers  are  distinct  from  those  of  uncon¬ 
trolled  rivers  in  several  respects.  Because  of 
main  stem  and  tributary  storage,  flow  increases 
are  more  gradual  during  a  flood  than  for  a  com¬ 
parable  uncontrolled  river.  Sharp-fronted  river 
waves  cannot  form  in  an  impoundment,  and  if 
wave  formation  occurs  upstream  or  in  a  tribu¬ 
tary,  attenuation  is  rapid  after  reaching  a  pool. 
Without  an  abrupt  wave  front,  the  development 
of  sufficient  stress  on  the  ice  cover  to  cause 
breakup  requires  locally  high  water  levels.  Also, 
the  flow  in  mainstream  backwater  reaches  oc¬ 
curs  at  greatly  reduced  energy  gradients,  corre¬ 
sponding  to  large  stress  reductions  and  an  in¬ 
creased  ability  of  the  ice  cover  to  resist  and  slow 
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the  progress  of  an  ice  breakup.  Together  these 
factors  typically  prevent  the  formation  of  a  sin¬ 
gle  breaking  front,  resulting  in  the  development 
of  breakup  segments  in  reaches  between  dams. 

Theory  of  dynamic  ice  breakup 

Ferrick  et  al.  (1986b)  presented  a  theory  of  dy¬ 
namic  ice  breakup  that  is  the  basis  for  our  analy¬ 
sis.  Briefly,  they  proposed  that  the  river  flow 
causes  hydraulic  forces  on  the  ice  cover  that  are 
resisted  by  the  ice  strength  and  points  of  sup¬ 
port.  The  primary  hydraulic  forces  on  the  ice  re¬ 
sult  from  the  shear  stress  of  the  flowing  water 
and  the  streamwise  component  of  the  weight  of 
the  cover.  A  good  approximation  of  the  applied 
hydraulic  force  per  unit  length  of  ice  sheet  is 

fh  =  BSj  (  y  R  +  Yjtj )  (1) 

where  Sf  =  flow  energy  gradient 
B  =  river  width  (m) 

Y  and  y,  =  specific  weights  of  water  and  ice 
(N/m3),  respectively 
R  =  hydraulic  radius  of  the  channel  (m) 
fj  =  ice  thickness  (m). 

The  applied  force  varies  linearly  with  the  energy 
gradient  of  the  river.  In  low-  to-moderate  gradi¬ 
ent  rivers,  the  energy  gradient  can  increase  by 
100%  or  more  on  the  front  of  a  river  wave  rela¬ 
tive  to  the  natural  stream  gradient.  In  steep,  shal¬ 
low  rivers  like  the  White  River,  the  increasing 
gradient  at  the  front  of  the  wave  is  accompanied 
by  large  increases  in  the  hydraulic  radius,  and 
the  combined  effect  produces  a  large  increase  in 
the  force  on  the  ice  cover.  A  positive  feedback 
exists  between  river  waves  and  ice  breakup.  Ice 
breakup  generates  unsteady  flow,  and  unsteady 
flow  increases  the  hydraulic  forces  on  the  sta¬ 
tionary  ice  cover,  potentially  leading  to  addition¬ 
al  breakup.  A  basic  requirement  of  an  analysis  of 
dynamic  breakup  is  the  capability  to  quantify 
this  interaction. 

An  analysis  of  ice  breakup  also  requires  the 
development  of  criteria  that  indicate  the  initia¬ 
tion  of  breakup.  Bolotnikov  (1982)  developed  a 
condition  based  on  the  theory  of  elasticity  for  the 
bending  failure  of  an  ice  cover  due  to  river  wave 
movement.  The  equations  given  for  both  the 
flexural  strength  and  the  breakup  resistance  of 
the  ice  cover  vary  linearly  with  ice  thickness, 
and  flexural  strength  is  the  critical  parameter 
that  determines  breakup  resistance.  We  have  es¬ 
timated  the  hydraulic  forces  at  failure  of  river  ice 


covers  and  have  found  them  to  result  in  distrib¬ 
uted  stresses  more  than  an  order  of  magnitude 
below  the  commonly  used  measures  of  ice 
strength.  Similarly,  Beltaos  (1985)  calculated  hy¬ 
draulic  forces  for  a  Thames  River  breakup  that 
were  also  much  less  than  the  usual  measures  of 
ice  strength.  We  conclude  that  dynamic  breakup 
occurs  either  as  a  result  of  a  failure  of  the  sup¬ 
port  of  the  cover  or  due  to  locally  high  forces 
that  result  from  ice  interaction  and  exceed  ice 
strength. 

The  ice  sheet  transfers  the  applied  forces  to 
the  banks  and  through  points  of  ice  grounding 
to  the  stream  bed.  The  support  provided  to  a  riv¬ 
er  ice  cover  by  the  bed  generally  decreases  with 
increasing  flow  depth,  and  the  primary'  supports 
are  at  the  banks.  The  hinge  cracks  present  at  the 
banks  are  continuous,  and  the  forces  on  the  ice 
must  be  transferred  across  these  cracks.  Prior  to 
breakup  the  applied  forces  are  in  equilibrium 
with  the  reactions  at  every  location.  The  primary 
resistance  to  ice  motion  frequently  relies  on 
crack  irregularity  and  the  interlocking  along  the 
crack  that  results.  When  two  surfaces  touch  at  a 
relatively  small  number  of  locations,  the  resis¬ 
tance  to  motion  is  proportional  to  the  actual  area 
of  contact,  which  varies  with  location  and  is  not 
generally  known.  It  would  be  reasonable  to  as¬ 
sume  that  the  contact  area,  like  the  surface  area 
of  the  crack,  increases  linearly  with  ice  thickness. 
Local  failures  of  the  ice  at  contact  points  are  re¬ 
quired  to  initiate  motion.  However,  the  stresses 
required  to  cause  the  failure  of  the  support  and 
the  motion  of  the  ice  cover  depend  on  several 
parameters  and  are  currently  unknown. 

At  locations  where  the  geometry  of  the  river 
does  not  permit  movement  along  the  hinge 
cracks  or  when  the  hydraulic  forces  are  in¬ 
adequate  to  cause  support  failure,  a  different  dy¬ 
namic  breakup  behavior  is  observed.  In  these  ice 
strength-dominated  cases,  locally  high  ice  forces 
that  overcome  the  flexural  strength  of  the  ice  are 
necessary  for  a  breakup  to  proceed.  The  under¬ 
turning  of  floating  ice  blocks  is  one  of  the  pro¬ 
cesses  of  interaction  between  brash  and  an  intact 
cover  that  produces  these  high  forces.  Daly  (in 
prep.)  developed  a  moment  balance  equation 
and  found  that  for  any  given  ice  thickness  and 
depth  the  tendency  for  block  underturning  in¬ 
creases  with  the  Froude  number.  Local  ice  fail¬ 
ure  would  be  expected  in  the  part  of  the  cross 
section  with  the  highest  Froude  number.  Other 
ice  interaction  processes  are  expected  to  respond 
similarly  to  flow  conditions. 
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The  energy  gradient  and  Froude  number  are 
closely  related  parameters.  Both  parameters  vary 
with  time  during  the  passage  of  a  river  wave, 
and  using  Chezv's  equation  the  relationship  can 
be  expressed  as 


v- 


c.\>x 


F2 

(C:/k) 
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where  F  =  Fronde  number 

C.  =  dimensionless  conveyance 
coefficient 

l7  and  V  =  average  velocity  (m/s)  and  depth 
(m)  in  a  cross  section,  respectively 
y  =  acceleration  due  to  gravity  (m/s2); 
k  relates  the  depth  and  hydraulic 
radius  of  wide  channels 
k  =  1  for  open  water  conditions  or 
k  =  2  if  the  channel  is  ice  covered. 

Ferrick  et  al.  (1986a,b)  studied  a  controlled  ice 
strength-dominated  breakup  and  noted  that 
downstream  progress  occurred  only  during  pe¬ 
riods  of  high  energy  gradient.  Therefore,  we  con¬ 
clude  that  a  parameterization  of  this  breakup  be¬ 
havior  could  be  expressed  as  a  function  of  either 
the  energy  gradient  or  the  Froude  number. 

Dynamic  ice  breakup  model 

A  one-dimensional  unsteady-flow  hydraulic 
model  is  the  basis  of  the  ice  breakup  model.  The 
open-channel  flow  continuity  and  momentum 
equations  are  solved  using  a  Preissmann  or  four- 
point  implicit  finite-difference  method  (Cunge  et 
al.  1980).  We  select  a  spatial  and  a  temporal  reso¬ 
lution,  supply  a  geometric  description  of  the  riv¬ 
er  (bed  slope,  channel  width,  ice  thickness),  and 
calibrate  the  hydraulic  roughness.  A  stationary 
ice  cover  increases  the  resistance  to  flow  of  the 
channei  relative  to  open  water.  When  an  ice  cov¬ 
er  is  present,  the  channel  roughness  parameter 
represents  a  composite  of  the  bed  and  the  ice, 
flow'  depth  is  measured  from  the  bed  to  the  un¬ 
derside  of  the  ice,  and  the  cross-sectional  area  of 
the  channel  is  determined  using  this  depth.  If  ice 
breakup  occurs,  (he  ice  is  assumed  to  move  at 
the  water  velocity,  and  flow  resistance  condi¬ 
tions  return  to  those  of  open  water. 

In  the  model,  breakup  is  evaluated  in  discrete 
reaches  that  are  set  by  the  spatial  resolution.  The 
ratio  of  the  spatial  to  temporal  resolution  is  the 


grid  celerity  c  ,  the  minimum  speed  of  the  break¬ 
up  in  the  model.  Because  of  the  feedback  be¬ 
tween  flow  dynamics  and  ice  breakup,  this  ratio 
must  be  less  than  or  equal  to  the  breaking  front 
speed  observed  in  the  river,  c^.  Otherwise,  an  ar¬ 
tificially  high  rate  of  water  release  from  storage, 
as  the  ice  breaks  upstream,  would  cause  a  more 
extensive  model  breakup  downstream  than 
would  occur  in  the  river.  The  speed  of  the  break¬ 
ing  front  can  be  used  to  define  a  breakup  Cour- 
ant  number  C  as 


r  -£t 'AL  -  £h 

A  A  o 


(3) 


I  he  breaking  front  speed  is  limited  bv  the  dy¬ 
namic  wave  celerity  c  =  V  +  \y Y  at  the  wave 
front.  The  dynamic  wave  celerity  appears  in  the 
usual  definition  of  the  Courant  number  Cr, 


(4) 


Finally,  we  obtain  a  pair  of  conditions  from  the 
relationships  between  c,„  q,  and  c  that  govern  the 
Courant  numbers, 


1  <  C  <  Cr.  (5) 

An  implicit  model  is  required  in  order  to  obtain 
both  a  stable  numerical  solution  and  the  correct 
rate  of  water  release  from  storage  when  the 
Courant  numbers  exceed  1.  In  the  case  of  a 
strength-dominated  breakup  where  q,  «  c,  eq  5 
requires  the  use  of  a  large  Cr. 

Available  theoiy  and  data  are  not  adequate  to 
quantitatively  model  the  failure  of  an  ice  cover 
during  strength-dominated  dynamic  breakup 
conditions.  However,  data  from  the  Connecticut 
River  (Ferrick  et  al.  1988)  are  adequate  for  an  ini¬ 
tial  evaluation  of  support-dominated  breakup  on 
that  river.  We  compute  the  hydraulic  forces  per 
unit  length  of  ice  cover  fh  at  each  model  section 
using  eq  1.  A  force  balance  written  for  the  ice 
cover  in  the  longitudinal  direction  yields 

Fn  dL  =  df  (6) 

where  Fn  =  fh  -  2 xtir  the  difference  between  the 
applied  forces  and  the  bank  reac¬ 
tions  per  unit  length  (N/m) 
dL  =  length  increment  in  the  longitudi¬ 
nal  direction  (m) 
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df  =  change  in  the  force  transmitted 

downstream  through  the  ice  cover 
corresponding  to  the  distance  dL 
(N) 

x  =  stress  transmitted  to  the  bank 
across  each  shore  crack  (Pa) 
xnl  =  maximum  allowable  stress  corre¬ 
sponding  to  failure  of  the  ice  cover 
support  (Pa),  x  <  xm. 

The  calculation  of  Fn  in  eq  6  assumes  equal  reac¬ 
tions  at  both  banks  to  balance  fh,  and  assumes 
that  these  forces  are  applied  over  an  area  of  a 
unit  length  times  the  ice  thickness.  An  ice  cover 
is  stable  when  the  forces  and  reactions  are  in  bal¬ 
ance  and  Fn  =  0.  As  the  hydraulic  forces  on  the 
ice  increase,  x  reaches  an  upper  limit  and  Fn  be¬ 
comes  positive,  indicating  increasing  forces  with 
distance  downstream.  Continued  ice  cover  sta¬ 
bility  requires  additional  support  to  counteract 
the  local  force  imbalance.  Ice  breakup  at  a  model 
section  occurs  when  x  exceeds  xm,  an  empirically 
determined  failure  condition  that  characterizes 
the  breakup  resistance  of  a  given  river  reach. 
The  capability  to  model  ice  breakup  in  regulated 
rivers  presents  several  options  for  ice  manage¬ 
ment,  including  control  of  breakup  by  flow  regu¬ 
lation.  The  remainder  of  this  discussion  devel¬ 
ops  a  case  study  of  controlled  breakup  on  the 
Connecticut  River. 

Control  of  Connecticut  River  ice  breakup 

The  flow  of  the  Connecticut  River  in  our 
study  reach  (Fig.  3)  is  controlled  by  Wilder  Dam 
upstream  and  Bellows  Falls  Dam  downstream. 
The  Connecticut  River  is  free-flowing  down¬ 
stream  of  Wilder  Dam,  with  an  average  bed 
slope  of  0.00037.  At  an  average  discharge  of  200 
nrVs  the  Connecticut  River  varies  between  100 
and  200  m  in  width  and  has  a  mean  depth  range 
between  1.5  and  3.0  m.  The  uncontrolled  White 
River  is  the  primary  tributary  in  the  reach,  enter¬ 
ing  about  2  km  downstream  of  the  dam.  The  lo¬ 
cation  of  the  head  of  the  Bellows  Falls  pool  var¬ 
ies  with  the  headwater  elevation  at  the  dam  and 
the  ice  conditions,  but  it  is  generally  near  data 
station  2  at  Windsor.  As  Wilder  Dam  does  not 
generally  pass  large  quantities  o'.  ..  ■,  the  White 
River  is  the  only  significant  i>  .  s.  'ce  at  break¬ 
up  external  to  the  study  reach  Because  of 

large  hourly  fluctuations  in  the  flow  release  at 
Wilder  Dam,  the  14-km  reach  betw  .  >  the  dam 
and  station  1  is  largely  free  >f  stable  ice  cover. 
However,  a  short  sinuous  reach  near  the  conflu¬ 


ence  with  the  Mascoma  River  generally  develops 
a  stable  ice  cover  that  is  resistant  to  breakup.  A 
continuous  stable  ice  cover  occurs  in  the  54-km 
reach  between  station  1  and  Bellows  Falls  Dam. 

In  extreme  dynamic  breakup  events  the  White 
River  rises  abruptly  to  a  high  peak  and  deposits 
large  quantities  of  ice  in  the  Connecticut.  Mean¬ 
while,  the  ice  on  the  Connecticut  River  is  compe¬ 
tent  and  intact,  and  the  combined  discharge  con¬ 
tinues  to  increase  rapidly  toward  a  peak  daily 
average  flow  in  excess  of  1200  nv'/s.  Under  these 
conditions  the  river  is  out  of  control  and  the 
probability  of  bridge  damage  or  loss  and  flood¬ 
ing  at  Windsor  is  high.  The  river  regulation  con¬ 
cept  applied  to  the  Connecticut  River  involves 
an  abrupt  flow  release  from  Wilder  Dam,  pat¬ 
terned  after  the  breakup  behavior  of  unregulated 
rivers,  with  a  minimum  initial  flow  and  a  mini¬ 
mum  Bellows  Falls  pool  elevation.  If  the  release 
has  a  sufficient  peak  discharge  and  duration,  the 
ice  cover  in  the  reach  upstream  of  station  3  will 
break  up.  These  controlled  conditions,  put  into 
effect  davs  in  advance  ot  the  White  River  ice  re¬ 
lease  and  uncontrolled  flows,  ensure  that  mini¬ 
mum  volumes  of  ice  and  \wiier  are  involved  in 
the  breakup,  minimizing  the  potential  for  ice 
and  related  flood  damage.  The  open  water  creat¬ 
ed  bv  the  late  winter  or  early  spring  breakup 
then  becomes  a  heat  source  that  collects  and  de¬ 
livers  heat  to  rapidly  melt  the  ice  accumulation 
downstream  of  station  2.  This  method  of  con¬ 
trolled  ice  breakup  would  not  produce  ice  break¬ 
ing  forces  farther  dowmstream  because  of  wave 
attenuation  in  the  Bellows  Falls  pool. 

The  Connecticut  River  regulation  required  to 
control  ice  breakup  in  this  reach  includes  the 
93.3-km  reach  between  Mclndoe  Falls  Dam  up¬ 
stream  and  Wilder  Dam,  and  the  68.2-km  reach 
between  Wilder  Dam  and  Bellows  Falls  Dam. 
Regulation  of  the  river  above  Wilder  Dam  is  nec¬ 
essary  because  the  available  storage  in  the  reser¬ 
voir  is  inadequate  to  produce  a  controlled  break¬ 
up  of  the  study  reach.  A  temporal  resolution  of 
0.5  hr  was  used  in  all  simulations,  and  the  uni¬ 
form  spatial  resolutions  used  were  3220  m  be¬ 
tween  Mclndoe  Falls  and  Wilder  dams,  and  2440 
m  between  Wilder  and  Bellows  Falls  dams.  The 
relatively  fine  spatial  resolution  below  Wilder 
Dam  yields  a  grid  celerity  of  half  the  observed 
speed  of  breakup.  Calibration  of  hydraulic 
roughness  of  the  Mclndoe  Falls-Wilder  segment 
of  the  model  was  achieved  by  comparison  with 
steadv-flow  water  surface  profile  data.  The 
Wilder-Bellows  Falls  model  segment  was  cali- 
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Figure  3.  This  reach  of  the  Connecticut  River  is  the  focus  of  the  ice  breakup  study, 
although  a  significantly  longer  reach  is  considered  in  the  model  simulations.  Data 
collection  stations  1,2,  and  3  are  indicated,  and  the  kilometer  points  measure  dis¬ 
tance  upstream  of  Bellows  Falls  Dam.  The  location  of  the  primary  damages  at  ice 
breakup  is  near  station  2. 
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brated  with  data  from  the  open-water  and  ice- 
covered  unsteady  flow  test  conditions  reported 
by  Ferrick  et  al.  (1988). 

Development  of  empirical  breakup  criteria 

The  rapidly  moving  support-dominated  dy¬ 
namic  ice  breakup  depicted  in  Figure  1  accurate¬ 
ly  represents  the  behavior  of  the  Connecticut 
River  in  the  study  reach.  Because  a  theory  that 
describes  the  support  failure  of  a  river  ice  cover 
is  not  available,  we  sought  empirical  failure  cri¬ 
teria.  The  approach  taken  was  to  apply  the  mod¬ 
el  to  simulate  the  dynamic  ice  breakup  of  27  Jan¬ 
uary  1986  reported  by  Ferrick  et  al.  (1988).  The 
river  flow  conditions  were  known  at  the  dams 
and  for  the  primary  tributaries,  and  the  ice 
breakup  behavior  was  observed.  Ice  thickness 
data  prior  to  breakup  were  not  available,  and 
uniform  thickness  was  assumed.  The  thicknesses 
of  ice  blocks  observed  at  station  2  were  about  0.3 
m.  Because  the  breakup  occurred  in  midwinter, 
the  ice  was  extremely  hard  and  characterized  by 
instantaneous  100-m-long  fractures  resulting 
from  collisions  between  large  ice  floes  and  a 
blunt  bridge  pier. 

The  failure  condition  at  each  section  in  the 
model  was  adjusted,  and  the  modeled  breakup 
wave  amplitude  and  subsidence,  breaking  front 
progression  and  speed,  and  ice  jam  location 
were  compared  with  observations.  The  basic 
character  of  the  observed  breakup  was  repro¬ 
duced  by  the  model,  with  breakup  stresses  on 
each  hinge  crack  exceeding  1.68  kPa  at  most  lo¬ 
cations  and  1.87  kPa  in  more  resistant  reaches. 


An  increased  resistance  to  breakup  may  have 
been  caused  by  greater  ice  thickness  in  reaches 
typically  subjected  to  larger  hydraulic  forces. 
These  stress  criteria  are  a  measure  of  ice  resis¬ 
tance  to  breakup.  We  assigned  the  stress  criteria 
of  this  initial  simulation  an  ice  resistance  of  1.0. 
The  results  of  this  simulation  compared  with 
field  data  in  Table  1  indicate  good  agreement  ex¬ 
cept  for  the  breakup  celerity.  Pariset  et  al.  (1966) 
developed  a  dimensionless  stability  diagram  for 
granular  river  ice  covers.  Their  results,  guided 
significantly  by  data,  indicate  that  the  value  of 
the  dimensionless  parameter  X  =  (B/YA)S,  must 
be  less  than  2.8x10'-'  for  an  ice  cover  to  resist 
breakup  at  any  ice  thickness  to  depth  ratio.  Con¬ 
sistent  with  intuition,  the  1.0  resistance  case  has 
X  =  7.4xl0'3  at  station  2  immediately  prior  to 
breakup,  representing  a  significantly  greater  re¬ 
sistance  to  movement  than  exhibited  by  granular 
ice  covers. 

These  celerity  and  stability  data  imply  that 
the  specified  1.0  breakup  stress  criteria  may  be 
somewhat  high.  Therefore,  the  breakup  stress 
and  ice  thickness  were  varied  for  the  same  flow 
conditions  to  both  improve  the  breakup  celerity 
agreement  and  determine  the  sensitivity  of  the 
simulation  to  these  parameters.  The  0.9  and  0.8 
ice  resistance  values  listed  in  Table  1  indicate 
10%  and  207c  uniform  reductions  in  the  stress  re¬ 
quired  at  breakup.  A  reduction  of  either  ice  resis¬ 
tance  or  ice  thickness  causes  an  increased  break¬ 
up  celerity.  Of  those  attempted,  the  0.9 
resistance,  0.30-m  thickness  case  yields  the  best 
agreement  with  the  available  data.  The  base  flow 


Table  1.  Comparison  of  observed  and  simulated  characteristics  of  Connecticut  River  ice  breakup 
near  station  2  on  27  January  1S36.  Wave  amplitude  and  subsidence  are  based  on  river  stage  1  hr 
prior  to  and  1  hr  after  the  peak  at  breakup,  respectively.  The  mean  energy  gradient  Sf  and  the  mean 
stress  at  breakup  were  determined  for  model  sections  between  stations  1  and  2.  Ice  jam  location  is 
the  distance  above  Bellows  Falls  Dam. 


Case 

(ice  resistance, 
ice  thickness) 

<m) 

Wave 

amplitude 

tint 

Wave 

subsidence 

(in) 

Average 
celerity  of 
breakup 
(m/s) 

Mean  S, 
lx  JO1) 

Mean  breakup 
stress 

(kPa) 

Ice  jam 
location 
l  km ) 

Breakup 

character 

Field  data 

0.3 

0.3 

2.7 

40.0 

Progressive 

1.0,  0.30 

0.33 

0.26 

1.35 

0.399 

1.80 

39.0 

Progressive 

1.0,  0.23 

0.16 

0.52 

0.362 

1.98 

17.0 

Simultaneous 

0.9,  0.30 

0.27 

0.19 

2.71 

0.377 

1.64 

39.0 

Progressive 

0.  8, 0.30 

0.16 

0.09 

0.362 

1.54 

26.8 

Simultaneous 
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in  the  river  exceeded  610  nrVs  in  the  study  reach 
prior  to  breakup.  The  modeled  peak  breakup 
discharge  at  station  2  averaged  860  rn'/s  with 
only  minor  variations  between  the  simulations,  a 
significant  increase  in  discharge  considering  the 
small  increase  in  river  stage  (Table  1).  The  break¬ 
up  in  the  reach  upstream  of  km  39.0  of  the  0.30- 
m  ice  thickness  occurred  at  an  increased  average 
stress  with  increasing  ice  resistance  conditions. 
In  this  reach  the  mean  energy  gradient  at  break¬ 
up  also  increased  with  ice  resistance,  marginally 
exceeding  the  mean  stream  bed  gradient  of 
0.00036  by  up  to  11%.  Without  high  gradients,  eq 
1  indicates  that  large  flow  depths  and  high  river 
stage  must  occur  to  produce  the  high  forces 
needed  for  a  dynamic  breakup,  and  very  high 
water  was  observed. 
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Figure  4.  Example  release  liydrograplis  of  Connecti¬ 
cut  River  dams  for  a  controlled  ice  breakup  experi¬ 
ment. 


Design  of  controlled  ice  breakup  experiment 

The  design  of  the  Connecticut  River  ice  break¬ 
up  experiment  had  to  satisfy  several  constraints. 
The  regulated  flows  must  comply  with  environ¬ 
mental  regulations  on  the  operation  of  the  river 
system.  Pool  fluctuations  required  by  the  test  are 
confined  to  the  ranges  of  normal  operation,  and 
minimum  flow  releases  are  maintained.  The  loss 
of  hydroelectric  power  production  and  other  ef¬ 
fects  on  normal  river  operations  are  minimized  if 
the  water  requirements  of  the  breakup  are  small. 
Together,  these  considerations  require  that  the 
hydraulic  forces  are  sufficient  to  cause  breakup 
at  a  significantly  lower  river  stage  and  smaller 
total  water  volume  than  occur  in  major  natural 
events.  The  flood  control  and  ice  damage  control 
objectives  are  ensured  by  the  small  volume  of 
the  release  and  by  excluding  the  ice  upstream 
and  in  the  White  River  from  participating  in  the 
breakup.  Excluding  the  ice  contributions  from 
the  tributaries  places  a  constraint  on  the  timing 
of  the  regulation.  The  need  to  prevent  additional 
ice  formation  following  the  breakup  requires  a 
net  heat  flux  into  the  river,  and  places  an  addi¬ 
tional  constraint  on  timing. 

Our  design  of  a  controlled  ice  breakup  on  a 
regulated  river  simulates  the  behavior  of  an  un¬ 
regulated  river.  An  abrupt  and  relatively  large 
short-duration  increase  in  flow  provides  a  con¬ 
trasting  condition  to  a  low-flow  initial  state  and 
produces  large  forces  on  the  ice  cover.  Proposed 
ice  breakup  flow  releases  at  the  three  Connecti¬ 
cut  River  dams  are  given  in  Figure  4;  Wilder  is 
the  only  dam  with  scheduled  releases  that  ex¬ 
ceed  turbine  capacity.  The  water  level  in  the 
Wilder  impoundment  must  initially  be  near  the 


top  of  its  operating  range  to  provide  the  capacity 
to  supply  the  release.  The  release  from  Mclndoe 
Falls  Dam  is  timed  with  the  Wilder  release  to 
maintain  the  Wilder  pool  in  its  normal  operating 
range,  slow  the  drawdown  near  the  dam,  and 
keep  the  segment  of  the  impoundment  experi¬ 
encing  rapid  drawdown  to  a  minimum.  The  Bel¬ 
lows  Falls  headwater  elevation  is  set  initially  at 
the  bottom  of  its  operating  range  to  move  the 
head  of  the  impoundment  downstream,  increas¬ 
ing  the  forces  attainable  downstream  of  station  2 
and  the  extent  of  the  breakup.  The  turbine  capac¬ 
ity  at  Bellows  Falls  Dam  is  sufficient  to  maintain 
that  elevation  with  only  minor  variations.  Stor¬ 
age  dams  on  a  pair  of  tributaries  above  (Om- 
pompanoosuc  River)  and  below  (Ottauquechee 
River)  Wilder  Dam  are  operated  to  minimize 
drawdown  and  enhance  the  breakup,  respective¬ 
ly- 

The  controlled  breakup  of  Connecticut  River 
ice  of  any  thickness  and  resistance  is  possible 
with  a  sufficient  release  volume.  However,  the 
usable  volume  severely  limits  the  release  and 
causes  an  uncertain  outcome  that  is  dependent 
on  the  properties  of  the  ice  cover.  A  sequence  of 
design  simulations  were  run  for  the  Figure  4 
flow  releases  with  a  range  of  ice  thickness  and 
resistance  parameters.  In  all  cases,  based  on  ob¬ 
servation  in  past  winters,  ice  covers  were  as¬ 
sumed  to  be  present  in  a  2.4-km  reach  below  the 
White  River  confluence  and  from  station  1  con¬ 
tinuously  to  Bellows  Falls  Dam.  A  summary  of 
important  characteristics  of  these  simulations  is 
presented  in  Table  2.  As  a  measure  of  conserva¬ 
tism  the  breakup  ice  stress  of  the  resistant  sec¬ 
tions  in  the  1.0  resistance  case  was  increased  to 
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Table  2.  Design  simulations  of  Connecticut  River  ice  breakup.  Wave  amplitude  and  peak  discharge 
are  at  station  2.  Average  breakup  celerity,  energy  gradient  Sf,  and  breakup  stress  were  determined 
for  model  sections  between  stations  1  and  2.  Ice  accumulation  location  is  the  distance  above  Bel¬ 
lows  Falls  Dam. 


CrtM* 

( ice  resistance, 
ice  thickness ) 

(tn) 

Pciik 

■ii<i  Inline 
(i  n  '/<) 

WjUV 

amplitude 

(m) 

Average 
celerity  of 
breakup 
(or  peak  stress) 
(mis) 

Alcan 
breakup 
or  peak 
s 

ixM  ) 

Mean 
breakup 
or  peak 
stress 

(Pn> 

hr 

ihXiiimiLUh’ii 

fiVii/ii'i; 

i  bin 

/>#'(  'akttp 
character 

1. 0,0.21 

725 

1.84 

2.26 

0.713 

2.69 

26.8 

Progress!  v e 

1.0,0.30 

761 

2.01 

1 .80 

0.676 

2.11 

34.1 

Progressive 

1.0,0.40 

501 

1.97 

(2.26) 

0.577 

1.50 

Non-breakup 

0.9,0.30 

746 

2.04 

1 .80 

0.717 

2.1  1 

26.8 

Progressive 

0.9,0.34 

746 

2.04 

1.80 

0.717 

1.92 

34.1 

Progress*  \e 

0.9,0.36 

741 

1.97 

0.600 

1.81 

3 1 .7 

Simultaneous 

0.9,0.40 

501 

1.97 

(2.26) 

0.577 

1.50 

Non -breakup 

0.8,0.40 

746 

2.04 

1.80 

0.717 

1.73 

34.1 

Progressive 

0.8,0.46 

501 

1.97 

(2.26) 

0.577 

1 .35 

- 

Non-breakup 

2.01  kPa.  The  standard  ice  stress  criterion  of  1.68  cases.  The  average  celeritv  of  the  other  progres- 

kPa  for  this  case  w’as  retained  from  the  January  sive  breakups  was  1.8  m/s  in  each  case.  The  si- 

1986  results.  As  before,  the  0.9  and  0.8  resistance  multaneous  breakup  was  the  most  resistant  case 

parameters  correspond  to  multipliers  that  re-  in  which  breakup  occurred,  and  represents  a 

duce  the  breakup  ice  stress  at  each  section  in  the  borderline  condition  between  the  breakup  and 

model.  non-breakup  regimes.  This  case  exhibited  late 

The  simulation  results  indicated  that  ice  and  concurrent  breakup  at  several  locations.  The 

breakup  would  occur  only  between  Wilder  Dam  average  of  the  peak  energy  gradients  at  model 

and  the  location  given  in  Table  2.  An  ice  breakup  sections  between  stations  1  and  2  was  much  larg- 

above  Wilder  was  not  predicted  in  any  of  these  er  in  all  of  these  design  simulations  than  in  the 

cases.  The  computed  peak  flow  velocity,  depth  January  1986  breakup  (Tables  1  and  2).  These 

and  discharge  at  all  locations  in  the  reach  were  high  gradients  provide  the  high  stresses  on  the 

less  than  or  equal  to  those  that  occurred  in  the  ice  sheet  at  the  relatively  low  river  stages  re- 

January  1986  breakup,  a  relatively  small  event  in  quired  in  a  controlled  breakup.  The  progressive 

the  historical  record  (Ferrick  et  al.  1988).  The  breakup  gradients  were  significantly  larger  and 

wave  amplitude  at  station  2  was  nearly  the  same  the  simultaneous  breakup  gradients  marginally 

for  all  cases.  However,  there  was  a  large  differ-  larger  than  the  non-breakup  gradients, 

ence  in  peak  discharge  between  the  progressive  The  flow  data  are  identical  for  the  cases  that 

breakup  and  non-breakup  cases.  The  resistance  did  not  generate  an  ice  breakup.  The  breakup 

of  the  ice  cover  to  breakup  for  a  given  case  is  cases  exhibit  minor  differences  in  stage-time  and 

characterized  by  the  product  of  the  resistance  discharge-time  response,  with  the  simultaneous 

parameter  and  the  ice  thickness,  corresponding  breakup  case  the  most  distinctive  of  the  group, 

to  an  applied  hydraulic  force  per  unit  length.  The  progressive  breakup  and  non-breakup  cas- 

The  0.21-m  ice  thickness  case  was  the  least  resist-  es,  compared  in  Figure  5,  reveal  several  impor- 

ant  to  breakup,  with  an  average  breakup  celerity  tant  differences  that  result  from  the  interaction 

between  stations  1  and  2  that  was  equal  to  the  of  the  ice  breakup  and  the  river  wave.  Without 

celerity  of  the  peak  stress  in  the  non-breakup  ice  breakup,  the  wave  amplitude,  the  peak  dis- 
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Figure  5.  Comparison  of  computed  stage  and  discharge  during  ice  breakup  and  non-breakup  simulations  at 
four  locations  below  Wilder  Dam. 


Figure  6.  Comparison  of  hinge  crack  stresses  on  a  0.4-m- 
thick  ice  cover  at  stations  2  and  3.  All  conditions  of  these 
two  cases  are  identical  except  for  a  reduction  in  the  resis¬ 
tance  of  the  supports  that  allowed  the  breakup.  In  this  case 
the  breakup  progressed  downstream  beyond  station  2  before 
stalling  about  a  kilometer  upstream  of  station  3. 


charge  and  the  rates  of  stage  and  discharge  in¬ 
crease  on  the  wave  front  all  diminish  with  dis¬ 
tance  downstream,  demonstrating  significant  at¬ 
tenuation  of  a  large  and  abrupt  river  wave  in  20 
to  30  km.  In  contrast,  ice  breakup  increases  the 
wave  amplitude  at  stations  2  and  3,  and  creates 
more  peaked  hydrographs  that  have  a  higher  or 
only  marginally  reduced  maximum  discharge 
relative  to  upstream  locations.  The  stage  and  dis¬ 
charge  increases  on  the  wave  front  have  a  much 
higher  rate  and  occur  earlier  at  these  stations  as 
a  result  of  ice  breakup.  All  of  these  tendencies 
emphasize  the  importance  of  the  release  of  water 
from  channel  storage  during  breakup  to  counter¬ 
act  wave  attenuation. 

Hinge  crack  stresses  at  stations  2  and  3  are 
compared  in  Figure  6  for  a  breakup  case  and  a 
non-breakup  case.  All  conditions  of  these  cases 
are  identical  except  for  a  reduction  in  the  resis¬ 
tance  of  the  ice  cover  supports  that  allowed  the 
breakup.  At  the  lower  resistance  the  ice  cover  at 
station  2  breaks  up,  but  the  breakup  stalls  just 
over  a  kilometer  upstream  of  station  3.  In  both 
cases  the  peak  stress  at  each  station  occurred 
about  0.5  hr  prior  to  the  peak  discharge  and 
about  1  hr  prior  to  the  peak  stage.  The  breakup 
causes  the  stresses  at  both  stations  to  increase 
earlier  and  attain  significantly,  higher  peak  val¬ 
ues  than  without  breakup.  The  short  duration  of 
the  high  stress  condition  caused  by  the  short- 
duration  release  and  the  rapid  attenuation  of  the 
peak  stress  immediately  downstream  of  the  stall 
prevent  the  formation  of  a  significant  ice  jam. 

Conclusions 

River  ice  breakup  is  complex,  exhibiting  wide 
variations  in  behavior  from  a  thermal  breakup  at 
low  forces,  resulting  from  a  large  decrease  in  ice 
strength,  to  a  dynamic  breakup  of  thick  and 
competent  ice  at  very  high  forces.  The  mode  of 
failure  of  an  ice  cover  in  a  dynamic  breakup  is 
determined  by  geometric  and  hydraulic  con¬ 
ditions,  varying  between  strength-dominated 
and  support-dominated  failure  at  different  loca¬ 
tions  in  a  river.  Support-failure,  strength-failure, 
and  ice  jamming  represent  a  continuum  of  dy¬ 
namic  ice  breakup  behavior.  These  processes 
typically  occur  in  sequence,  giving  dynamic 
breakup  a  surge-stall  character.  The  develop¬ 
ment  of  abrupt,  high-amplitude  river  waves  dur¬ 
ing  a  breakup  on  an  unregulated  river  requires 
significant  runoff  and  competent  ice  cover.  Reg¬ 
ulated  rivers  respond  differently  at  breakup  than 


unregulated  rivers  because  the  control  structures 
change  both  the  hydraulic  conditions  and  the  ice 
conditions. 

The  fundamental  component  of  the  theory 
presented  is  the  intrinsic  relationship  between 
unsteady  flow  and  ice  breakup  for  dynamic 
breakup  conditions  with  a  rapidly  moving 
breaking  front.  A  dynamic  ice  breakup  model 
that  treats  support-dominated  failure  was  devel¬ 
oped  and  applied  to  the  Connecticut  River.  Data 
from  a  dynamic  breakup  were  used  to  obtain 
empirical  criteria  for  the  failure  of  the  support  of 
the  ice  cover.  The  relationship  between  river 
waves  and  ice  breakup  was  observed  in  field 
data  from  the  White  River  and  was  demonstrat¬ 
ed  in  the  model  simulations  of  the  Connecticut 
River.  The  simulations  also  indicate  that  the 
breakup  response  of  a  river  changes  from  pro¬ 
gressive  to  simultaneous  and  then  to  non¬ 
breakup  conditions  for  the  same  initial  wave  by 
increasing  the  ice  thickness  or  breakup  resis¬ 
tance.  These  results  identify  the  need  for  an  in- 
situ  index  test  of  ice  cover  resistance.  Further  re¬ 
finement  and  verification  of  our  theory  of  ice 
breakup  would  be  achieved  most  readily  by  con¬ 
ducting  controlled  breakup  field  experiments.  A 
high  spatial  density  of  stage-time  records  and  ice 
observations  with  known  boundary  conditions 
are  needed  for  several  cases  to  adequately  char¬ 
acterize  both  strength-dominated  and  support- 
dominated  dynamic  breakups. 

Our  theory  of  dynamic  breakup  provides  the 
framework  for  implementing  controlled  breakup 
and  minimizing  potential  ice-related  damages 
on  regulated  rivers.  Control  of  ice  breakup  by 
river  regulation  uses  existing  dams,  requires  a 
relatively  small  volume  of  water,  and  is  environ¬ 
mentally  sound.  The  concept  is  to  remove  the  ice 
cover  from  upstream  of  locations  with  a  high 
damage  potential  during  breakup,  and  to  en¬ 
hance  the  rate  of  melting  of  the  ice.  The  rapid 
stage  increase  and  short-duration  characteristics 
of  a  controlled  breakup  are  patterned  after  those 
of  unregulated  rivers.  Because  of  the  short  dura¬ 
tion  of  the  release,  the  ice  forces  at  the  down¬ 
stream  extent  of  the  breakup  rapidly  diminish, 
preventing  the  formation  of  a  significant  ice  jam. 
The  involvement  of  a  limited  ice  volume  and  the 
low  initial  river  stage  and  discharge  distinguish 
this  event  from  uncontrolled  breakup  events  on 
regulated  rivers.  Once  developed  for  a  river  sys¬ 
tem,  controlled  ice  breakup  can  be  implemented 
on  relatively  short  notice  and  only  when  needed. 
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below. 


Ferrick,  Michael  G. 

Framework  for  control  of  dynamic  ice  breakup  by  river  regulation  /  by  Mi¬ 
chael  G.  Ferrick.  Hanover,  N.H.:  U.S.  Army  Cold  Regions  Research  and  Engi¬ 
neering  Laboratory;  Springfield,  Va.:  available  from  National  Technical  Informa¬ 
tion  Service,  1989. 

vi,  21  p.,  illus.,  28  cm.  (CRREL  Report  89-12.) 

Bibliography:  p.  14. 

1.  Controlled  ice  breakup.  2.  Dynamic  ice  breakup.  3.  Flood  control.  4.  Ice  jams. 
5.  Ice  management.  6.  River  ice.  7.  Unsteady  flow.  I.  Mulherin,  Nathan  D.  II.  Unit¬ 
ed  States  Army.  III.  Corps  of  Engineers.  IV.  Cold  Regions  Research  and  Engineer¬ 
ing  Laboratory.  V.  Series:  CRREL  Report  89-12. 


